Effects of Rescattering in (e, e'p) Reactions within a Semiclassical Model 
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The contribution of rescattering to final state interactions in (e, e'p) cross sections is studied 
for medium and high missing energies using a semiclassical model. This approach considers two- 
step processes that lead to the emission of both nucleons. The effects of nuclear transparency 
are accounted for in a Glauber inspired approach and the dispersion effects of the medium at low 
energies are included, ft is found that rescattering is strongly reduced in parallel kinematics. At 
high missing energies and momenta, the distortion of the short-range correlated tail of the spectral 
function is dominated by a rearrangement of that strength itself. In perpendicular kinematics, a 
further enhancement of the experimental yield is due to strength that is originally in the mean field 
region. This contribution becomes negligible at large missing momenta. 

PACS numbers: 25.30.Fj, 25.30.Dh, 21.60.-n, 21.10.Pc, 21.10.Jx. 



I. INTRODUCTION 



Nuclear correlations strongly influence the dynamics of 
nuclear systems 0, 0- In particular, the repulsive core 
at small internucleon distances has the effect of depleting 
the shell model orbitals and inducing high momentum 
components in the nuclear wave functions 3, 4]. The 
main effects of short-range correlations (SRC) consist in 
shifting a sizable amount of spectral strength, about 10- 
15% j5j, to very high missing energies and momenta, to- 
gether with increasing the binding energy |^ . The result- 
ing reduction of the occupation numbers of the deeply 
bound orbitals appears to be fairly independent of the 
given subshell and of the size of the nucleus, except for 
a slight increase with the central density of the system. 
Theoretical studies of the distribution of short-range cor- 
related nucleons for finite nuclei have been carried out 
using a local density approximation (LDA) by Benhar et 
al. I3 and with many-body Green's functions by Miither 
et al. 0. These calculations suggest that most of the 
missing strength is found along a ridge in the momentum- 
energy plane {pm-Em) which spans several hundreds of 
MeV/c (and MeV). Such behavior is confirmed by recent 
experimental data f§|. 

It is important to note that the depletion of spectro- 
scopic factors for closed shell orbitals observed near the 
Fermi energy is more substantial than the 15% reduction 
discussed above due to long range effects like the cou- 
pling to collective modes We note that this reduc- 



* Electronic 



address: 



lbarbieri@triumf.cal 



tion also tends to be weaker for loosely bound orbitals 
like halo states [Tol [lll|. As a consequence, it becomes 
particularly interesting to study the spectral distribution 
in heavy nuclei where the mean field single particle or- 
bitals extend to regions far from the Fermi level and tend 
to decouple from surface effects. A measurement of the 
spectral function for the complete mean field region of 
^°^Ph has been undertaken recently at NIKHEF T^l. In 
this region the single particle orbitals are sensibly frag- 
mented and it is required to probe missing energies up 
to 100 MeV or more. Further detailed information 
on SRC could be obtained from (e, e'p) experiments that 
directly search for the missing strength at very high miss- 
ing momenta and energies. This is particularly appealing 
since the details of its distribution strongly infiuence the 
binding energy of finite nuclei and nuclear matter ^6^ . Be- 
sides, it could shed more light on how much the interior 
of large nuclei is sensitive to the effects of finite size and 
proton-neutron asymmetry. 

Unfortunately, past measurements of the short-range 
correlated tail by means of (e, e'p) reactions have been 
limited due to the enormous background that is gener- 
ated by final state interactions (FSI), see for example 
Refs. 0, The issue of how to minimize the FSI 

has been recently addressed in Ref. There, it was 

suggested that FSI in exclusive (e, e'p) cross sections are 
dominated by two-step rescattering processes like the one 
depicted in Fig. This becomes particularly relevant 
when regions of small spectral strength are probed in 
perpendicular kinematics |39j |. A study of several kine- 
matic conditions shows that the rescattered nucleons can 
move spectral strength in the Pm-Em plane, from the top 
of the ridge toward regions where the correlated strength 
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FIG. 1: (Color online) Schematic representation of the di- 
rect knockout of a proton (a), given by the PWIA, and the 
contribution from a two-step rescattering (b). In the latter 
a proton or neutron is emitted with momentum pa and dif- 
ferent missing energy and momentum (_E^,p^). Due to a 
successive collision, a proton is eventually detected with the 
same momentum p/ seen in the direct process. 



background noise. Other possible contributions that in- 
volve the excitation of a A resonance are expected to be 
more sensitive to transverse degrees of freedom. Parallel 
kinematics tend to be more clean due to the high momen- 
tum that is required for the detected proton. In Ref. (l^ , 
it was suggested that the contribution of rescattering can 
be diminished by choosing parallel kinematics and tak- 
ing advantage of modern electron beam facilities. New 
data were subsequently taken in these conditions by the 
E97-006 collaboration at Jefferson Lab H E [13 for a 
set of nuclei ranging from carbon to gold. Clearly, FSI 
could still play a role even in parallel kinematics and need 
to be properly addressed before the relevant physical in- 
formation is extracted from the experiment. We note 
that a similar dependency of the FSI on the kinematics 
is also predicted in Ref. ^3 for (e, e'NN) reactions in 
superparallel kinematics. 

The issue of computing the effects of rescattering 
has been considered recently by means of the multi- 
step dynamics approach and by using Glauber the- 
ory j23,|2ll|23|. All these calculations suggest that multi- 
ple rescattering contributions (more than two-steps) are 
relatively small in light nuclei like ^^C but can become 
relevant for large systems. Interference effects between 
FSI and SRC correlations can also play a role I2H . How- 
ever, all these effects were seen to be reduced in parallel 
kinematics. In Ref. ^^Ii the scattered proton was de- 
tected at energies at which a full distorted wave calcula- 



tion, in terms of an optical potential, is required. How- 
ever, rescattering processes leading to the emission of two 
nucleons (one of which is not detected) can lead to the 
reappearance of a part of the experimental strength ab- 
sorbed by inelastic processes. This effect was investigated 
in terms of a semiclassical model inspired by the work of 
Ref. [23. Even if very different kinematical situations 
were considered, the reaction mechanism included in the 
latter approach is the same as pointed out in Ref. 
and Fig. ^ Some partial simplifications occur for high 
energy protons, since the relevant effects of the medium 
are limited to Pauli blocking. Therefore this approach 
offers a valid starting point to investigate the FSI effects 
needed in the analysis of data at large missing energy 
and momenta. Other effects such as meson exchange 
currents and the excitation of resonances also need to 
be investigated. However, these are beyond the scope of 
the present paper and will be considered in future work. 
In this paper, we consider the approach of Ref. jl^ and 
extend it to high missing energies. We then apply it to 
the kinematics of both the NIKHEF T^l ad the E97- 
006 [T^ experiments to evaluate the importance of 
two-step processes for the different kinematics employed. 

The model for computing the contribution of rescat- 
tering is depicted in Sec.^ together with a discussion of 
the inclusion of the absorption effects in terms of nuclear 
transparency. A practical application requires the knowl- 
edge of the in-medium differential cross section, which 
is calculated in Sec. IIIII by extending the approach of 
Ref. Sees. WJ\ and HVBl report on the results for 

the kinematics used in the above experiments, at medium 
and high missing energies respectively. Conclusions are 
drawn in Sec. IVl 



II. MODEL 

This work considers contributions to the experimental 
yield that come from two-step mechanisms in which a 
reaction (e, e'a) is followed by a scattering process from 
a nuclcon in the medium, N'{a,p)N" , eventually leading 
to the emission of the detected proton. In general, a 
may represent a nucleon or another possible intermediate 
particle. In this work we will only consider the channels 
in which a is either a proton (with N' = p or n) or a 
neutron. In the following we will also use the letter a to 
label the possible open channels. 

Following the semiclassical approach of Refs. 0, , 
the contribution to the cross section coming from rescat- 
tering through the above channels is written as 
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dEa dQi. dEr dQp 



xPT(Pa;ri,r2)p7V'(r2) PT(p/;r2,oo) , (1) 



where {Eo^'ko) and {Ef,pf) represent the four- moment a 
of the detected electron and proton, respectively. Eq. ^ 
assumes that the intermediate particle a is generated in 
plane wave impulse approximation (PWIA) by the elec- 
tromagnetic current at a point ri inside the nucleus. Here 
-ftT = |pa|-Ea is a phase space factor, S^{p'j^, E!^)/M is the 
spectral function of the hit particle a normalized to one 
[i.e., M = N{Z) if a is a neutron (proton)] and a^^^ the 
off shell electron-nucleon cross section, for which we have 
used the ccl prescription of de Forest . The pair dis- 
tribution functions gaN'Qi^i — r2|) account for the joint 
probability of finding a nucleon N' in r2 after the parti- 
cle a has been struck at ri |125| . The integration over the 
kinetic energy Ta of the intermediate particle a ranges 
from to the energy u transfered by the electron. The 
nuclear transparency factor Pt{p', Ti, r2) gives the trans- 
mission probability that the struck particle a propagates, 
without any interactions, to a second point r2, where it 
scatters from the nucleon N' with cross section (JaN' ■ 
The whole process is depicted in Fig. ^p. The point nu- 
cleon densities PAr(r) are normalized to either the num- 
ber of neutrons or protons. These were derived from ex- 
perimental charge distributions by unfolding the proton 
size '26]. We employed equal distributions for neutrons 
and protons, which is a sufficiently accurate approxima- 
tion since even in ^'^'^Pb neutron and proton radii differ 
by less than 4%. 

The nuclear transparency for the ejected nucleon was 
considered in Ref. _23] according to Glauber theory. The 
probability Pt that a proton struck at ri will travel with 
momentum p to the point r2 without being rescattered 
is given by 

Pt(p; ri,r2) = 
exp{- dz [5pp(|ri - r|) CTpp(p,p(r)) pp(r) (2) 



(|ri - r|) CTp„(p,p(r)) p„(r)] } , 

where the z axis is chosen along the direction of prop- 
agation p, an impact parameter b is defined so that 
r — b-|-zp, and zi (^2) refer to the initial (final) position. 
Eq. ^ differs from the standard Glauber theory by the 
inclusion of the pair distribution functions {\ri — r|). 
In principle, the gpN functions should depend on the den- 
sity and on the direction of the inter-particle distance. 
However, these effects has been shown to be negligible 
in Ref. |23l |. In the present application we find that a 
simple two-gaussian parameterization of the ^pjv can ad- 



at saturation density. The in medium total cross sec- 
tions (Tpp {p, p) and dpn (jp, p) used in this work have been 
computed according to Ref. ^"^^ account for the ef- 
fects of Pauli blocking, Fermi spreading and the effective 
mass generated by the nuclear mean field. For energies 
above 300 MeV they have been extended to incorporate 
the effects of pion emission . 

The nuclear transparency is defined as the average over 
the nucleus of the probability that the struck proton 
emerges from the nucleus without any collision. This 
is related to Pt by 



r= dTpp{y)PT{p;v,^) 



(3) 



It should be mentioned that the practical experimental 
definition of nuclear transparency depends on the spe- 
cific kinematics employed and that Eqs. (0) and Q rig- 
orously apply only to the parallel case 28]. For the 
present approach, these are the right quantities to be 
included in Eq. ^ since they describe the loss of flux in 
the direction of propagatiom. In the case of ^°*Pb and 
an outgoing proton with energy Ef ~ 1.1 GeV (kinetic 
energy of -161 MeV), Eq. © gives T = 0.37. With 
Ef — 1.8 GeV, which is of interest for the calculations 
of Sec UVBl T = 0.63 for and T = 0.29 for ^^^Au. 
The PWIA contribution of the direct process. Fig. 
also needs to be corrected for the nuclear transparency 
effects, 



d^ 



0'PlV/A(T) 



dE^ dilf, dEf dflp 



KaTp' S';{p'^,E'JT , (4) 



with T given by Eq. (PJ. Eq. Q is consistent with the 
assumptions of Eq. from which it would be obtained 
if the rescattering event was substituted with the limit of 
r2 to infinity. 



III. EVALUATION OF THE IN-MEDIUM 
NUCLEON-NUCLEON RATE 

The last ingredient required in Eq. is the in medium 
cross section for the rescattering process, in which the 
particle a (either a proton or a neutron) hits against a 
bound nucleon in its way out, eventually leading to the 
emission of the detected proton. This can be evaluated 
by extending the approach of Ref. to describe the 
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the NN differential elastic cross section in free space 



dcr. 



pN 



\fpN{cos 61)1' 



\M{s,t,u)\ 
64 7r2 s 



where s, t and u are the Mandelstam invariants {^/s being 

energy in the cm. system) and Al(s,i,u) represents 
the square of the Lorentz invariant amplitude 29] aver- 
aged (summed) over the initial (final) spins. For nucleon 
momenta above 1 GeV/c and for small angles, the scat- 
tering amplitude in Eq. jSJl is well approximated by its 
central part and can be written as 



fpN 



^ -riot 

An 



(epN + i)exp{~Plj^{p^ - Pff/2} , (6) 



where pi and p/ are the initial and final momentum of 
the scattered nucleon, EpN is the ratio of the real to imag- 
inary part of the scattering amplitude and cr*'jj is the 
total scattering cross section. At low energy the values 
of A4{s,t,u) where extracted from the SAID phase shift 
data analysis [s^l- For the pp case, we chose to keep the 
differential cross section constant for angles smaller than 
5° and larger than 175°, in order to avoid the Coulomb 
peak in the forward and backward directions. However, 
the results of the present work are largely insensitive to 
this choice of the cut off angle due to Pauli exclusion. The 
solutions of the SAID program were used for energies in 
the laboratory system up to 1.6 GeV for pp scattering 
and 1.2 GeV for the pn case, which are well contained 
within the range of validity of this data base. At higher 
energies Eq. © was used with parameters (7*°* =44.0 mb. 



=41.1 mb and Cj, 



-0.48. The slope coeffi- 



cients were chosen by requiring that Eq. (O yields the 
correct ' 
implies 



correct values ct^'^ for the total elastic cross section. This 



pN 



2 N tot 2 
''pN) "pN 



16 7T af^ 



where and tT^Jj were extracted from the experi- 
ment l3lf . A direct comparison for energies above 1 GeV 
showed that Eq. © appropriately approximates the data 

(5) 

from the SAID data base, thus there exists an overlap re- 
gion where these approaches are both accurate and join 
smoothly. 



In deriving the rescattering rate we assume that the in- 
teraction between the two nucleons is localized enough so 
that the amplitude M[s, t, u) is not altered by the pres- 
ence of the surrounding nucleons. However, the medium 
can sensibly modify the cross section due to the spectral 
distribution of the momentum of the hit nucleon and due 
to the effects of Pauli blocking. In the spirit of the local 
density approximation (LDA), which underlies Eq. 
the momentum of the hit nucleon, iV', is taken to be lo- 
cally distributed as in infinite nuclear matter. The den- 
sity of the latter being the one of the point Y2 where the 
collision occurs, pnm = Pn'{^2)- The effects of the nu- 
clear surface are eventually included by integrating over 
PN'{'^2)dY2 in Eq. QJ. For the present purposes it is 
appropriate to further approximate the symmetric nu- 
clear matter with a free Fermi gas |40j . The assumption 
of a completely filled Fermi sea is also consistent with 
the Dirac-Brueckner-Hartree-Fock (DBHF) employed be- 
low. Initially, the hit nucleon N' is in the Fermi sea 
and therefore must have a momentum h smaller than 
kp = [iir'^pissM/'^Y^^- At the same time the Pauli prin- 
ciple requires that the particles in the final state will have 
momenta pf and 1, both larger than kp. Among all the 
nucleons involved in the process, p/ refers to the detected 
proton while the others can be either neutrons or protons 
depending on the channel a. The probability per unit 
time of an event leading to the emission of a proton with 
momentum p/ is obtained by imposing the Pauli con- 
straints and integrating over the unobserved momenta h 
(7) and 1, 



N' 



dpf d^pj 
= 2p}0{pf-k 



2 9{pf - kp) 



9{kF - h)e{i - kp) Wi 



{2nf 



(2^ 



:e{kp - h)0{l - kp) 



1 



\M{s, t, u) 



647r2 Ea{pa)EN'{h)Efipf)EN"il) 
I 



6{Ea+EN' -Ef-EN") 



l=Pa+h-p/ 



where is the volume of a normalization box and, for the pion production threshold is related to Eq. (|SJ) by 



a free nucleon, E^ip) = (p^ + m^)^/^^ jj^ 



is the probability per unit time for the event + ft.^ — > 
-I- , which can be expressed in terms of M.{s,t,u). 
The inverse life time of the nucleon a for energies below 



1 

Ta 



E 

N'—p,n 



dpf 



aN' 



dpf dQ 



(9) 



At low energies, a nucleon traveling through the 
medium acquires an effective mass due to dispersion ef- 
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scalar field Us and the time component of a vector field 
Uv |32| . This particular approach allows to maintain the 
relativistic framework adopted in Eq. (|SJ|. The values of 
Us and t/y in nuclear matter were computed in Ref. [s^ 
by solving the DBHF equations. The results were found 
to be consistent with the value of the non relativistic ef- 
fective mass extracted from nucleon-nucleus scattering. 
At the energies considered by the DBHF calculations, 
Us and Uv are predicted to be essentially momentum 
independent. Thus, the energy of a nucleon moving with 
momentum p is given by 



mDipNAi) = riiN + Us{pnm) 



(10) 

En{p,Pnm) = \Jp'^ + rajjipNAi) + Uv{pnm) ,(11) 

where mu is the Dirac effective mass. It should be noted 
that Us and Uy have large and opposite values. The suc- 
cess of the DBHF approach rely on a subtle cancellation 
of their effects and require a self-consistent calculation of 
the interaction with the medium |33j |. As a consequence, 
the use of Eas. HlU|) and to compute the in medium 
cross section at low energy does not guarantee a priori 
accurate predictions. In a relativistic model, this would 
require a more elaborate calculation of the scattering am- 
plitude M (see for example Refs. 0,|33). In the present 
work, we are simply interested in giving an approximate 
treatment of the dispersion effects on the rescattering 
of protons for the kinematics that will be considered 
in Sec. IIV Al The kinematics relevant for studying the 
short-range correlated tail involve much higher energies, 
where the nuclear cross section is known to approach the 
free one. Above 1 GeV, Us and Uy arc not known except 
that they are no longer momentum independent and that 
they should decrease to zero. In this case, Pauli blocking 
gives the only relevant contribution of the medium and 
the dispersion effects are negligible. Therefore Us and 
Uv will be set to zero in the calculations of Sec. IIVBI 

When the effective mass is accounted for, the momen- 
tum of a nucleon participating in the rescattering pro- 
cess is related to its energy E{p) by Eq. ifTTl) . In general, 
Al(s, t, u) is off shell. However, following the assumption 
that the interaction is not appreciably modified by the 
in-medium effects we use the on shell values extracted 
from the vacuum pN cross section, Eq. In vacuum, 
this depends on only two invariants that were chosen to 
be s = (p^ + /i^)^ and t — {p'^ ~ and computed 
accounting for the dispersion relation l|ll|) . The energy 
denominators appearing in Eq. ((SJ were also taken to be 
equal to the total energy of the nucleon, Eq. (fTT|) . We 
note that this differs from the normalization of a Dirac 
spinor in the medium js^ . However, this prescription is 
consistent with the choice of using the free scattering am- 
plitude M.{s,t^u) since it provides the right normaliza- 
tion in the non relativistic limit. Finally, the in medium 
scattering rate is given by 
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FIG. 2: (Color online) In medium elastic pN cross sections 
as a function of density for laboratory energies of 161 MeV 
(top panel) and 250 MeV (bottom panel). The curves show 
the vacuum cross section (full line), the results obtained by 
accounting for Fermi spreading and Pauli blocking (dashed 
line) and by adding the interaction with the medium through 
Eq. mill (dot-dashed line) . 



where p^i is the density of the hit nucleon and Vg — 

dEa{pa)/dpa = Val \/ vi + "f^^Da the group vclocity of 
the incoming one. Note that the dispersion effects mod- 
ify the rescattering rate, Eq. p2l) . in three different ways. 
First, both Vg and the Jacobian dpf/dEf depend on Us 
and Uv- Second, the density of final states for the scat- 
tered nucleons is modified by using Eq. (|lll) in the energy 
delta function 5{Ea{pa) + E^'iK) - E}{pf) - En„{1)). 
Third, the energies of nucleons just above kp are low- 
ered by Eq. Hll|) , which allows scattering at energies that 
would otherwise be Pauli forbidden. 

Figure |2 shows the in medium effective cross section 
as a function of density obtained by integrating Eq. I|12|) 
over angle and energy. Two values of the energy of the 
incoming nucleon are considered. The solid line gives the 
vacuum cross section while the dashed line includes the 
effects of Pauli blocking and Fermi spreading. A further 
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FIG. 3: (Color online) In medium elastic pN cross sections 
computed for normal nuclear matter density as a function of 
the nucleon energy. The curves show the vacuum cross section 
(full line), the results obtained by including Pauli blocking 
and Fermi spreading (dashed line) and also the dispersion 
effects (dot-dashed line). 



fects. While the calculation of the Pauli blocking effects 
is equivalent to the work of Ref . i23i] , the full cross section 
obtained here is somewhat smaller than the one obtained 
in the corresponding non relativistic result. A difference 
between the two approaches should be expected since the 
present calculation is based on the Dirac effective mass, 
given by Eq. (|10|l . This is different form the non relativis- 
tic definition of effective mass, which is instead related to 
the vector potential Uy 36j . The effects of Pauli blocking 
at higher energies can be seen in Fig. O where the elastic 
cross section is computed at normal nuclear matter den- 
sity pnm = 0.16 fm~'^ for energies up to 3 GeV. As it 
can be seen, the effects of Pauli blocking remain relevant 
at large energies where they produce a constant reduc- 
tion of the cross section. The further reduction due to 
the effective mass affects the results at low energies and 
tends to become less important at 1 GeV even for con- 
stant values oi Us and Uy- However, it is meaningless to 
extend the calculation of the effective mass effects above 
this energy since the values of Us and Uy are unknown 
in this region. 



IV. RESULTS 

Eq. requires the knowledge of the undistorted spec- 
tral function S^{pm, Em) of the target nucleus. For the 
present purposes, the strength in the mean field region 
can be described as a sum over the nuclear orbitals 
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FIG. 4: (Color online) Theoretical results for the rescatter- 
ing contribution to the reduced spectral strength of ^°*Pb 
for the kinematics of Ref. p^. The full (dashed) lines re- 
fer to the kinematics with lower (-Bo ) and higher (E^) en- 
ergy beams. The black lines shows the input spectral func- 
tion, Eq. 11311 . The results obtained for rescattering with two 
nucleons emitted in the continuum are given by the green 
curves. For the kinematics E^, the plots show the effects 
due only to Pauli blocking and Fermi spreading (dot-dashed 
lines). If the in medium dispersion was accounted for but no 
two nucleon emission was imposed, Eq. Q would give to the 
dot-dot-dashed curves. 
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(13) 



where $i(pm) are the single particle wave functions, Zi 
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variable width |37 



fi{En 



V{Er, 



a {E„ 



Ep 



(14) 



(15) 



b+iE^- EfY ' 

where the Fermi energy was taken to be Ep — 6 MeV, 
a = 24 MeV and b = 500 MeV^. The neutron spectral 
function was also obtained from Eq. (|13l) by including all 
the neutron orbits occupied in the shell structure. For Pb 
and Au no experimental data are available for neutrons, 
thus the occupation numbers Zi were assigned by extrap- 
olating the trend measured for protons in ^^"Pb ^^l- For 
kinematics chosen to probe medium and low missing en- 
ergies, as those discussed in Sec. lIV Al S'l,jp{pm, Em) cov- 
ers all experimentally accessible strength. In Sec. IIVBI 
this will need to be extended by including the distribu- 
tion of nucleons in the SRC correlated region. 

In the following, the rescattering yield computed from 
Eq. Q has been converted to a reduced spectral function 
by dividing it by \pfEf\Taff^ , evaluated for the kinemat- 
ics of the direct process [according to Eq. Q)]. This gives 
the straightforward correction to the model spectral func- 
tion that is due only to single rescattering effects. 



A. Application to ^°®Pb 

The studies of Ref. ; employed two different parallel 
kinematics in which the outgoing proton was emitted in 
the same direction as the momentum transfer (thus q 
and Pm were also parallel). The central kinetic energy 
of the outgoing proton was kept constant at 161 MeV 
{pf =570 MeV/c) in both cases. The two kinematics 
differ only for the energy of the electron beam, which 
was taken to be E^ —674 MeV for the first case and 
EJ^ =461 MeV for the other. These choices correspond 
to a virtuality ranging between 0.08 and 0.22 GeV^. 

In applying the model of Sees. and IIIII to (e, e'p) 
reactions one has to impose the further constraint that 
both nucleons are emitted in the continuum. This is triv- 
ially satisfied for high energy nucleons, for which the in- 
teraction with the medium can be neglected. When the 
dispersion effects are included from Eq. (|11() one has to 
require that Ei{l) > in the integrand of Eq. ||SJ. Pro- 
cesses in which the undetected nucleon remains bound 
are beyond the scope of the present work and would re- 
quire a proper quantum mechanical treatment to include 
its reabsorption vertex. At low energies this should be 
analyzed in terms of a proper optical potential model. 

The results of Eq. ^ for the rescattering contributions 
leading to two nucleons in the continuum are plotted in 
Fig. ^ The yield resulting from rescattering is between 
one and two orders of magnitude smaller than the direct 
signal, except for low missing momenta and missing ener- 
gies above 60 MeV, where it gives a correction of about 
20%. The rescattering effects are also found to be in- 



For comparison we also show the results obtained by in- 
cluding only Pauli blocking and Fermi spreading and the 
ones obtained when the in medium effects are included 
without requiring that the undetected nucleon is in the 
continuum [i.e. allowing Ei{l) < in Eq. JSJ]. Although 
the latter result has no direct physical relevance, it shows 
that only a part of the suppression of the rescattering ef- 
fects is due to the corrections for effective mass in the 
medium ('to_d' curves of Fig. The remaining reduc- 
tion is a consequence of the energy required to reach the 
two-nucleon emission threshold. 



B. Proton knock out from the SRC region 

This section considers the results for the kinematics of 
Ref. 0, where the aim is to directly probe SRC. In this 
case it is convenient to write the spectral function as the 
sum of a mean field and a correlated part, 

S'^ iPrn, Em) = S'ijp{pm,Em) + StrriVm, E^) , (16) 

where S^gj.j.{pm, E^,) describes the short-range correlated 
tail at very high missing energies and momenta 0, Q ■ In 
the present work this was parametrized as 



^corr ^Pm 5 En 



C e-"P" 



[Em - e{pm)? + [r(p™)/2]2 



(17) 



where e{pm) and T{j)m) are smooth functions of the miss- 
ing momentum that were chosen to give an appropriate 
fit to the available ^^C(e,e'p) data in parallel kinemat- 
ics The solid line in Fig. [3 shows the model spec- 
tral function, Eq. (|16|l . employed in the present calcula- 
tions for that part of the Pm-Em plane where S'^^^r domi- 
nates. The calculation with a ^^^Au target employed the 
same ^^^^ of Eq. lO multiphed by 79/6 or 118/6 to 
account for the correct number of protons and neutrons, 
respectively. This is shown for protons in Fig. |H1 At 
energies close to the Fermi level the hole spectral func- 
tion is dominated by its mean field component. For ^^C 
these are orbitals in the s and p shells, which are known 
experimentally and represent about 60 % of the total 
strength. The spectroscopic factors and wave functions 
used in Eqs. H13fl -(|15ll are the ones extracted from the 
world data in Ref. [38|. Since no direct data are avail- 
able for gold we choose to employ the spectral function 
discussed above for the neighbor nucleus ^°^Pb 12] but 
modifying the occupation of the last orbitals to account 
for the different number of nucleons in those shells. 

We have performed calculations of the rescattering 
contributions by employing same sets of three parallel 
and two perpendicular kinematics used in Ref. |16 |. In 
the parallel case, the initial momentum of the proton 
(~Pm) was centered at different angles with respect to 



the momentum transfered by the electron, dq 



25, 



21 and 36 deg, while the corresponding energies of the 
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FIG. 5: (Color online) Theoretical results for the total (direct 
+ rescattered) reduced spectral strength in the correlated re- 
gion. The results are given for parallel (dashed line) and per- 
pendicular kinematics (dot-dashed line). The full lines show 
the model spectral function, Eq. (I16II . employed in the calcu- 
lations. All panels refer to a ^'^C target and employ the same 
line convention. Note that the results for different sets of par- 
allel kinematics do not always overlap exactly. This is mostly 
due to the dependence of the off-shell cross section a^"}} on 
the kinematics [l^ . 



This implied angles -dqf ^ 0.5, 1.5 and 10 deg, respec- 
tively, between the final proton and the momentum trans- 
fered. For the two perpendicular kinematics, the same 
angle -dqpi ^ 90 deg was used, while Ef ^ 1.25 and 
1.35 GeV and -dgf ^ 25 and 29 deg. The four momen- 
tum transfered by the electron was always in the range 
~ 0.3-0.4 GeV^. More details on these kinematics 
are discussed in Ref. [l^ . 

Due to the loss of energy of the ejected nucleon at 
the rescattering vertex, the spectral strength is always 
shifted toward higher missing energies. This is clearly 
visible in the results for both ^^C and ^^^Au, which are 
shown in Figs. Eland |Sl(the sum of direct plus rescatter- 
ing signals is plotted). The contribution to parallel kine- 
matics is negligible at missing energies below the peak of 
the correlated tail but it tends to become more impor- 
tant for Em > 150-200 MeV. This confirms the expected 
trend that a part of the strength seen in this region is 
dragged from places where the hole spectral function is 
larger [T5l |. The same behavior is seen in perpendicu- 
lar kinematics where, however, rescattering effects are 
already relevant at small missing energies. In this situa- 
tion the direct process accounts for only 30-50% of the 
total yield obtained at the top of the correlated peak. 




FIG. 6; (Color online) Theoretical results for the total re- 
duced spectral strength in the correlated region. The results 
are given for parallel (dashed line) and perpendicular kine- 
matics (dot-dashed line) . The full lines show the model spec- 
tral function, Eq. 11611 . employed in the calculations. All pan- 
els refer to a ^^'^Au target and employ the same line conven- 
tion. Note that the results for different sets of kinematics do 
not always overlap exactly. This is mostly due to the depen- 
dence of the off-shell cross section on the kinematics [iql . 



PWIA signal by more than an order of magnitude. It 
should be noted that for both parallel and perpendicular 
kinematics the FSI become more important as the mass 
number increases. In general, this is due to the average 
distance that the outgoing nucleon has to travel inside 
the nucleus. Thus it carries a dependence on the nuclear 
radius. Figure [7| compares the results for both target 
nuclei in parallel kinematics. For comparison the yield 
for ^^^Au has been normalized to the same number of 
protons of ^^C. The reduced spectral function of gold is 
indeed always larger. The same consideration applies for 
perpendicular kinematics in Fig. [S] 

To study the origin of the rescattered strength, the cal- 
culations were repeated for gold by neglecting the mean 
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FIG. 7: (Color online) Theoretical results for the total re- 
duced spectral strength of ^^C obtained in parallel kinemat- 
ics (dashed line) compared to the analogous results for ^^^Au 
(dot-dashed line), normalized to the number of protons of 
carbon. The full line shows the input spectral function of 
Eq. 11611 employed in the calculations. All panels employ the 
same line convention. 



mean field strength extends to large missing energies, up 
to ~100 MeV. Thus rescattering effects can easily spread 
it into the correlated region. Figure IHl compares the the- 
oretical reduced spectral strength of Au with the anal- 
ogous result obtained when S'ljp is included. As one 
can see, relevant contributions from the mean field ap- 
pear for momenta up to about 500 MeV/c. The results 
at higher missing momenta are completely dominated by 
the rescattering from the correlated tail S^g„ into the 
correlated region itself. The situation is instead differ- 
ent in parallel kinematics where the present results for 
two-step rescattering do not shift any strength from the 
mean field region, even for a large nucleus like Au. Such 
large shift appears to be energetically forbidden do to 
the large energy of the scattered proton adopted in these 
kinematics. 



V. CONCLUSIONS 

A proper understanding of the relevance of FSI and 
their dependence on the kinematics is important in mod- 
ern (e, e'p) experiments that attempt to observe the cor- 
related strength at medium and high missing energies. 
The present work suggests a semiclassical approach to 
compute the effects of two-step rescattering, which is one 
of the leading contributions at high proton energies, and 
applies it to investigate its consequences for the kinemat- 




FIG. 8; (Color online) Theoretical results for the total re- 
duced spectral strength of ^'^C obtained in perpendicular kine- 
matics (dashed line) compared to the analogous results for 
^^'^Au (dot-dashed line), normalized to the number of protons 
of carbon. The full line shows the input spectral function of 
Eq. I|16|l employed in the calculations. All panels employ the 
same line convention. 



The model assumes a PWIA for the electromagnetic 
vertex, in which the struck nucleon is described by the 
full hole spectral distribution of the target nucleus. This 
gives the possibility of investigating how FSI shift the 
original strength from the direct process within the miss- 
ing energy and momentum plane. The absorption effects 
of the medium were accounted for by means of trans- 
parency factors. The rescattering is described in terms of 
the differential nucleon-nucleon cross section modified in 
order to account for Pauli blocking and Fermi spreading 
effects. The dispersion effects due to the nuclear medium 
have been included from DBHF results at the energies 
where they are relevant. 

For kinematics involving outgoing protons of the order 
of few hundreds of MeV, the present model was employed 
to estimate the reappearance of strength through inelas- 
tic channels that lead to two-nucleon emission. In the 
reaction ^''^Pb(e, e'p) the overall effects were seen to be 
more than an order of magnitude smaller than the direct 
signal, slightly increasing for small missing momenta and 
missing energies above 60 MeV. This supports an anal- 
ysis of the experimental data based on usual distorted 
wave calculations. 

The same model was applied for the kinematics of the 
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FIG. 9: (Color online) Reduced spectral strength of ^^'^Au 
computed in perpendicular kinematics as generated by the 
full spectral function of Eq. 1)16^ (dashed line) or by the sole 
correlated part S^orr (dot-dashed line). The full line shows 
the model spectral function of Eq. 1161 . Its mean field compo- 
nent, S%jp, is not visible in this plot except for small missing 
momenta. Note that the dashed and dot-dashed lines over- 
lap in the bottom panel. All panels employ the same line 
convention. 



distributions at high momenta. Calculations were per- 
formed for ^^C and ^^^Au targets, which have different 
radii. In general, rescattering was found to be much 



smaller in parallel kinematics than in perpendicular ones. 
In the latter case a large amount of strength is shifted 
from regions where the spectral function is big to regions 
where it is smaller, thus overwhelming the experimental 
yield from the direct process. This confirms the studies 
of Ref . [31 ■ The contribution from rescattering effects is 
also seen to increase with the nuclear radius. The rescat- 
tering of nucleons originally emitted form the mean field 
orbitals was found to be important in perpendicular kine- 
matics and for missing momenta lower than ~500 MeV/c. 
No such a large shift of strength was found in the par- 
allel case. The remaining effects of rescattering are due 
to a rearrangement of the spectral strength within the 
correlated tail itself. 

The present results provide a good first insight in the 
redistribution of strength due to FSI in (e, e'p) reactions. 
However, it is clear that in order to properly explain the 
real experimental yield observed at high missing ener- 
gies and for heavy nuclei other effects beyond the two- 
step rescattering need to be addressed [23, El. Rele- 
vant extensions of the present formalism would include 
corrections from meson exchange and nucleon excitation 
processes, as well as an investigation of the importance 
of multiple rescattering for heavy nuclei in parallel kine- 
matics. Work in this direction is in progress. 
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